Population control of 2s-2p transitions in hydrogen 
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We consider the time evolution of the occupation probabilities for the 2s — 2p transition in a 
hydrogen atom interacting with an external field, V(t). A two-state model and a dipole approxima- 
tion are used. In the case of degenerate energy levels an analytical solution of the time-dependent 
Shrodinger equation for the probability amplitudes exists. The form of the solution allows one to 
choose the ratio of the field amplitude to its frequency that leads to temporal trapping of electrons 
in specific states. The analytic solution is valid when the separation of the energy levels is small 
compared to the energy of the interacting radiation. 

PACS numbers: 32.80.Qk,42.50.-p 



I. INTRODUCTION 

The ability to control quantum systems, e.g. to cause 
a complete transfer of population from an initial state 
to a desired final state usingan external field, is crucial 
in quantum computation an d modification of reac- 
tion dynamics Since the 1970s many techniques 
for controlling single quantum systems have been devel- 
oped, e.g. methods for trapping atoms ,.5] and arranging 
atoms on surfaces |6j. In this paper we consider transi- 
tions between the 2s and the 2p levels in hydrogen caused 
by an interaction with a microwave radiation field. Hy- 
drogen is an instructive example, since it is the simplest 
naturally occuring atom, and consequently one of the few 
that permit detailed and precise comparison with theo- 
retical models. If one limits the number of active states 
to two, one has a qubit. 

Transitions between two discrete states can be, under 
certain conditions, accurately described within the two- 
state model H H ES IH El 13. The energy sepa- 
ration between the 2s and 2p levels, discovered in 1947 
by Lamb and Rctherford [Tj], is small compared to the 
nearest levels, Is and 3p. This enables one to use the 
two-state model over a wide range of field frequencies. 
The interaction between a two-state atom and a single 
mode radiation field is often treated using the rotating 
wave approximation (RWA) that makes the problem ex- 
actly solvable 0, M, 13 El- It was first discussed in 
detail in 1963 by Jaynes and Cummings 01 ■ Although 
RWA provides an accurate description in many cases of 
interest, it has its limitations as well. It requires the 
frequency of the radiation field to be close to the transi- 
tion frequency, and it cannot be applied to a system with 
nearly degenerate energy levels. 

Another analytic solution of the time-dependent 
Shrodinger equation exists in the limit of degenerate en- 
ergy levels. We present this solution here. We show that 
it has an interesting feature which can be used in the 
problem of quantum control. 

Detailed analysis is given of the applicability of our 
analytical solution. Calculations are presented for 2s — 2p 
transitions in atomic hydrogen. 

Atomic units are used throughout the paper. 



II. FORMALISM 

In the two-state approximation the exact time- 
dependent wave function has the form 



¥(?,t) = ai(f)Si(f) + 02(t)*2(f) 



(1) 



where are the orbital functions, and a n are the 

probability amplitudes. The time variation of the prob- 
ability amplitudes is determined by the time-dependent 
Schrodinger equation, namely, 
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ia 2 
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V 2 iax 



(2) 



where V$j = / $*(r)V ext (r, t)$j(r) d 3 r, and E n are the 
electronic energy levels. The diagonal matrix elements 
Vu and V 22 are required to be zero by the dipole selection 
rule. 

In the dipole limit, an external electromagnetic field 
interacts with an electron via the time-dependent poten- 
tial 



V ext (f,t) = -rE{t) , 



(3) 



where r is a relative electron-nuclear distance, and E is 
an electric field which for the plane monochromatic wave 
can be written as E = IEq cos (uit). The matrix elements 
Vij can then be expressed as 

Vtj(t) = -fijcEo cos (wt) , (4) 

where f tJ — f* % = J §*(r)r$j(r) d 3 r. We define w 2 i = 
E 2 — Ei , and introduce the Rabi frequency 

X = r 21 eE . 



(5) 



The arbitrary zero of energy is set at E\ , so E 2 — > E 2 — 
E\ = u>2i- Then equations for the probability amplitudes 
become 



iai — — X cos(tut)a 2 

ia 2 = uj 2 \a 2 — xcos(ujt)ai 



(6) 



These equations can be solved numerically for an arbi- 
trary values of energy separation ui 2 i between electronic 
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energy levels, a field frequency u>, and a Rabi frequency 
X that expresses the field-atom interaction in frequency 
units. 

Analytical solutions exist in various limits. When both 
lu and u>2i are large compared to the frequency of oscil- 
lation of the transition amplitudes and are close one to 
another, an analytical solution can be obtained either 
under RWA or within the oscillating-field theory [lif . 
Equations © decouple when either % — > or u>2i — > 00. 
Replacing an oscillating field by a stationary field (i.e. 
taking the limit u> — ► 0) also leads to an analytically 
solvable problem. Another analytic solution, considered 
in detail below, exists in the limit of degenerate energy 
levels. 



A. Analytic solution: limit L021 — ► 

Equations (JSJ) can be solved analytically in the limit 
^21^0, namely 

a\ = cos[— sin(u;f)] 

ui 

a,2 = i sin[— sin(wi)] . (7) 

CO 

Here we take the initial conditions as ai(0) = 1 and 
fl 2(0) = , i.e. initially the system is in state 1. 

As we can see from (JTJl , the temporal behavior of the 
probability amplitudes in this limit is defined by the ratio 
x/cj. The argument of the "outer" function (e.g., sin for 
02) varies between ±x/k>. Thus, the transition amplitude 
oscillates between and sin(x/w). Note that if x/w < 
7r/2, the amplitude 0,2 (and, consequently, the occupation 
probability P%) never reaches 1, i.e. complete population 
of the second level cannot be achieved. If x/w > tt/2, 
the amplitudes oscillate in a quasi-regular fashion. 

B. Complete population transfer 

An interesting solution occurs when the ratio xl 1 ^ = 
7r/2. Then Eqs (0 take the form 

a\ = C0S [~^ sin(cjt)] 

a,2 = z sin[— sin(u;i)] . (8) 

The corresponding occupation probabilities Pi — \a\\ 2 
and P2 = |ft2| 2 are 

Pi = cos2 [7^ sin(wi)] 

P 2 = sin 2 [|sin(wi)] . (9) 

For this special value of xl 10 the probability is completely 
transfered back and forth between the two levels with 
the period tt/uj as shown in figure ^ It can be easily 
shown that when reaches its maximum value at 




0.2 0.4 0.6 0.8 1 

Time (in units of the field period, T=2ji/co) 

FIG. 1: Analytical solution for x/ 1 - = 7r /2 for a system with 
two degenerate energy levels. 

t = 7t/2cj, a small deviation e from the value 7r/2 in the 
ratio x/ w reduces P2 from 1 to 1 — e 2 . This can be used 
to estimate how sensitive this approach is to fluctuations 
in the strength of the external field, x- The period of 
oscillation of the occupation probabilities is two times 
smaller than the period of oscillation of the radiation 
field since Pi — \a.i\ 2 . 

Complete population transfer occurs in general wher- 
ever the action, A(t) = J* V2i(t) g?t, is an integer multi- 
ple of ft/4. Each cycle has a "flat" part where the prob- 
ability value remains very close to the extreme values of 
or 1 for an extended period of time (cf. figure 0. As 
we can see, the external radiation field with the ratio 
x/lo = ir/2 can be used to cause periodical population 
inversion of the electronic levels (a physical realization 
of a quantum controllable system). An analysis of the 
applicability of this analytic solution to 2s — 2p transi- 
tions in hydrogen, including the criteria for the choice of 
a field frequency, is given in Section ^QJ Modifying the 
shape of V(t) to further control population transfer is 
also discussed. 



C. Duration of the population transfer 

Equations (© can be expanded in a power series in the 
vicinity of the point to = T/4, where Pzit) reaches its first 
maximum. Then P2 can be accurately approximated by 
taking only first few terms in the series, which converges 
rapidly for r < T = 2ir/w. For the form of the potential 
chosen in Eq.(J3J), the first, second and third derivatives 
of the function P2 are zero at t = to, and 

P 2 (to + r) = l + I^V + 0(r 6 ) « . (10) 

Hence the occupation probability P 2 in the vicinity of its 
maximum can be approximated by a 4th degree polyno- 
mial. 
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FIG. 2: Population leakage in the vicinity of to = T/4. Solid 
line - analytic solution ©, dash-dotted line - polynomial ap- 
proximation il l lit . 

The polynomial approximation of Eq. (|10f) tells one how 
to choose the field frequency to provide a desired duration 
of the populated state. To obtain a duration T s of the 
populated state with a population leakage less than some 
critical value P cr , the field frequency to be used is 

Near r = the error is 0((lot) 6 ). This is illustrated in 
figure [3 



III. RESULTS 

In this section we compare full and analytic results 
for the 2s — 2p transition in hydrogen in the two state 
approximation. 



A. Population leakage 

To obtain our analytic solution to Eq.©, we had to 
neglect the ^2102 term in the equation ia,2 = ^210-2 — 
XCOs(u>t)ai. This approximation requires that U21 <C Xj 
or (since x/ w = 7r /2) I ^21 *C lo. When LO21 is finite, the 
population transfer is not complete. The effect of this 
population leakage can be calculated by expanding the 
transition amplitude 02 in a power series in time includ- 
ing terms u>2i- It is easily shown that the first 0*21 term 
corresponds to the second derivative, djj(t). The differ- 
ence between the exact and the analytic solutions for the 
transition amplitudes is Aa,2(t) ~ iu^iX^ 2 ? an d 

AP(t) = \Aa 2 (t)\ 2 « \uj\ iX ^ . (12) 
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FIG. 3: Occupation probability P2(t) in the vicinity of to = 
T/4: analytical solution, Eq.©, and numerical calculations 
for different values of the ratio lo/lovi (1, 10 and 100). De- 
viation from the analytic solution decreases as (1021 /lo) 2 (cf. 
Eqn. O l 

By the time the occupation probability P2 reaches its first 
maximum at to — T/4 = tt/2lu, the difference becomes 

AP{t Q ) ~\C-f£fY . (13) 

Comparison of Eqs. l|13fl an d (|11|) shows that the choice 
of the field frequency, lo, involves a trade-off between the 
duration of the state and the population leakage. The 
higher the frequency, lo, the smaller the population leak- 
age, AP. But the higher u>, the shorter the duration T s 
of the population transfer. 

B. Calculations of 2s 2p transitions in hydrogen 

Here we compare the full (numerical) calculations and 
the analytical solution for a 2s — 2p transition in hydro- 
gen. The full calculations were done by numerically in- 
tegrating Eqs. © using a standard fourth order Runge- 
Kutta algorithm. In hydrogen the energy separation of 
2s and 2p levels (Lamb shift) is 4.37* l0~ 6 eV, while the 
next available level (3p) is ss 1.89eV away. Therefore, 
one can choose the frequency of the external field lo such 
that LO21 <C lo <C u>2s3p, i-e. both the degenerate level 
approximation (limit LO21 — > 0) and the two-state model 
can be used. 

As shown in figure [21 the difference between the full 
and analytic solutions for Pi (t) = 1 — P2 (t) in the vicin- 
ity of to = tt/2lo is large when lo = 1021- But it de- 
creases for lo > L02i- The difference is less than 1% for 
L0/LO21 = 10, and 0.01% for lo/lo 2 i = 100. Since that dif- 
ference is proportional to (0021 /lo) 2 , one can go up in 
frequency as high as fractions of eV (that brings the 
difference between the approximate and exact solution 
down to 10~ 10 — 10~ n ) and yet be far sensibly from the 
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nearest available (2s — 3p) resonant transition frequency. 
Therefore, a radiation field with the wavelength from a 
few yitm (and the intensity of the order l0 12 W/cm 2 ) up 
to a few cm (and the intensity of the order 10 4 W/cm 2 ) 
can be used for 2s — 2p transitions in hydrogen. The 
probability of the multiphoton excitation to the 3p state 
(or any other state, including continuum) also seems to 
be small for the range of the external field frequencies 
under consideration. 

In summary, our analytic approximation appears to 
be valid for 2 s — 2p transitions in hydrogen over a broad 
range of field frequencies. 



C. Changing the shape of V(t) 

We have noted previously that the choice of the field 
frequency is a trade-off between two competing factors: 
reducing the population leakage and increasing the du- 
ration of the populated state. If wishes to obtain a long 
lasting populated state with very small leakage, it may 
be that both requirements for the system cannot be met 
simultaneously for a single frequency radiation field. In 
this case it may be possible to use another form of the 
interaction potential. In particular, one may change the 
shape of the "flat" part of the probability (cf. figure^) by 
using different shapes for the external potential V ex t(t) . 
For an arbitrary external potential the formulas for the 
occupation probabilities can be written in the form: 



Pi(t) = cos 2 [/ V21 (r)dr] 



P 2 (t) = sin 2 [ / V 21 {t)<1t] . 
Jo 



(14) 



A Taylor series expansion can be used for choosing the 
shape of the external potential. For example, one can 
use the potential with the first non-vanishing derivative 
of the order higher than four to make the shape of the 
populated state even flatter. The derivatives can be cal- 
culated using the general formula for the nth derivative 
of a composite function , which in this case takes the 
form 



d n P 2 
dt n 



n \ d m F y' . y" . . 



i\j\...k\ dy m v lr v 2! 



where P 2 {t) = F{y) = sin 2 y , y = y(t) = f Q V 21 (r) dr , 
^2 indicates summation over all solutions in non-negative 
integers of the equation i + 2j + • • • + Ik = n and m = 

i + j H h k. 

In principle, Ean. (|15fl tells us how to shape V(t) to 
control the population transfer. As one can see from 
Ea. (|15fl . one may eliminate all terms up to the order r k 
by choosing the potential for which all derivatives up to 
the order (k — 1) are zero at the point t = t . That will 
"flatten" the shape of the populated (or depopulated) 
state, i.e. allow one to use smaller frequency (and, as a 



result, an increased duration of the state) to achieve the 
same level of population leakage. Therefore, the shape of 
an external potential can be used (along with the choice 
of the field frequency) for quantum control. For example, 
in a truly two-level system, choosing the potential of the 
form 



V(t) = -S(t - t ) 



leads to 



Pi(t) 



l-e(i-to) 

0(t-*o) • 



(16) 



This represents complete and immediate population in- 
version at t = U). However, Fourier transformation of 
V{t) now contains all frequencies, so one can no longer 
use a two-state approximation for the 2s — 2p transition 
in hydrogen, since higher states become necessarily in- 
volved. 



IV. DISCUSSION 

In this paper we have used the 2s — 2p transition in 
hydrogen as an example of population transfer. The same 
approach can be used for any other atomic or molecular 
system that has a similar pattern of energy levels (two 
levels located close one to another and far from the other 
levels). As we have shown in Ea. HlHJI . for an external 
field with a single frequency w, the difference between 
the exact calculations and the analytic approximation 
varies as (cj 2 i/^) 2 for lo 2 i/lo -C 1, and does not depend 
on the internal structure of the atom or molecule. This 
feature opens the possibility for using different systems 
with different values of the transition frequency, uj 2 \ and 
different ranges for the field frequency, lo. 

A downside of analyzing more complex atomic or 
molecular systems is that analytic expressions for the 
orbital functions (and, consequently, for the matrix el- 
ement V21 in Eq. I|14[)) are not always available. Then 
genetic algorithms (GA) may be used to choose a shape 
for the external potential, V(t). The application of GA 
with active feedback to the selective breaking and mak- 
ing chemical bonds in polyatomic molecules has been dis- 
cussed in detail by Rabitz et al p^. |20| , and tested ex- 
perimentally [2lL l22|). As in any optimization scheme, 
the effectiveness of the GA increases significantly when 
the initial value of the parameter close to the optimal one 
is chosen. One might combine our method with GA us- 
ing approximate analytic orbital functions to choose the 
"starting" form of the potential, estimate T s and P cr , 
and then employ the GA scheme. 

For 2s — 2p transitions in hydrogen, sources of the 
microwave radiation with corresponding intensities may 
now be available, which can be used to test our model. 
Metastable excited H(2s) atoms could interact with a 
low frequency (e.g. microwave) radiation. Since the life- 
time of the 2p state is much shorter than the 2s state, 
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the population of the 2p state may be monitored by ob- 
serving photons emitted in 2p — Is transitions. In these 
experiments the smallness of the Lamb shift requires the 
use of the temperatures below the 1 mK to exclude the 
thermal transitions between the levels. 

Finally we note that the approximation introduced 
here is related to the removal of time ordering in the 
time evolution operator 23]. The full evolution opera- 
tor may be expressed as U(t, to) — Texp(— i J V(t')dt'), 
where T is the Dyson time ordering operator that puts 
the sequence V"(t n )...V(£2)V"(fi) in the order of increas- 
ing time. When time ordering is removed T — > 1 and 
U(t,to) — > exp(— i f t V(t')dt'). In a two state system 
this corresponds to a\ = cos(f V(t')dt') and a 2 = 

isin(f* V(t')dt'), as used in this paper. This approxi- 
mation has been recently studied in the context of sec- 
ond order perturbation theory 0, HH > where it has 
been shown that the time evolution of different electrons 
becomes uncorrelated as T — ► 1. Corrections to this 
Magnus-like approximation |27| lead to correlated prop- 
agation that is non local in time. These corrections cor- 
respond to quantum fluctuations in the energy in short 
lived intermediate states arising from coupling of a quan- 



tum system with its macroscopic environment |2; 



V. SUMMARY 

A new analytical solution of two coupled channel equa- 
tions has been found that enables one to temporally con- 
trol the electron population of 2s and 2p states in hydro- 
gen by using a time varying external field. In addition, 
the population leakage can both be easily estimated, and 
than be further controlled by changing the shape of the 
external field. Our method can be applied to various 
systems. 
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